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■ ABSTRACT 

On , The Sachs- Wolfe effect is known to produce large angular scale fluctuations in the Cosmic Microwave 

■ Background Radiation (CMBR) due to gravitational potential fluctuations. We show how the angular 
correlation function of the CMBR can be expressed explicitly in terms of the mass autocorrelation function 

' ^(r) in the Universe. We derive analytic expressions for the angular correlation function and its multipole 

^ , moments in terms of integrals over ^(r) or its second moment, Jsir), which does not need to satisfy the sort 

of integral constraint that ^(r) must. We derive similar expressions for bulk flow velocity in terms of ^ and 

(N ■ 

J3. One interesting result that emerges directly from this analysis is that, for angles 0, there is a substantial 
contribution to the correlation function from a wide range of distance r and that the radial shape of this 
I contribution does not vary greatly with angle. 

^ ' Subject Headings: Cosmic Microwave Background — Cosmology: Large Scale Structure of the Universe. 

(N : 

cn . 

' 1 Introduction. 

. One sensible explanation for the currently clumpy state of the Universe is that it arose from the grav- 

' itational amplification of an initially nonuniform background mass density. The anisotropy of the Cosmic 

. Microwave Background Radiation (CMBR) provides a signature of these primordial density perturbations. 

In simple terms, variations in the gravitational potential arising from density fluctuations at the last scatter- 
ing surface produce variations in the redshift of photons across the sky. This is the Sachs- Wolfe effect (Sachs 
I ' and Wolfe 1967, Wolfe 1968), which is the dominant contributor to CMBR fluctuations on large angular 

2 . scales. 

' The CMBR anisotropy is most easily quantified in terms of the angular correlation function, C{d) = 

Ci ■ {ST16T2) / (T)'^ , of the temperature averaged over all directions in the sky; here 6 is the angle between the 

^ ' directions 1 and 2. For potential fluctuations in a perturbed, flat, Friedmann universe, Sachs and Wolfe 

•i-H . (1967) related C{d) to the present power spectrum of matter fluctuations, P{k), by following photon paths 

' from the last scattering surface to the observer. Since the power spectrum itself is just given by the Fourier 

transform of the mass autocorrelation function ^(r), it is also possible to express C{6) in terms of the 
autocorrelation function. This was done for some special cases by Traschen and Eardley (1986). 

It is also well known that there are several other effects which can and do contribute to the CMBR 
anisotropy. These include Dopplcr shifts and thcrmalization effects associated with the surface of last 
scattering. However, on large angular scales, the dominant contributor to fluctuations is the Sachs- Wolfe 
effect. 

In discussing the angular correlations of the CMBR, it is traditional to subtract out the contribution to 
C{9) arising from its monopole and dipolc moments (cf. Martinez-Gonzalez and Sanz 1989; Gorski 1991). 
The monopole term is not physically interesting, while the dipole contribution to the observed correlation of 
temperature is usually assumed to be dominated by the peculiar motion of the observer with respect to the 
local standard of rest. In that case the dipole velocity is not itself directly related to potential fluctuations. 
The recent positive detection by COBE of fluctuations in the CMBR (Smoot et al. , 1992) has been reported 
with monopole, dipole and quadrupole moments of C{9) removed to exclude these large contributions to the 
anisotropy from sources other than the Sachs- Wolfe effect. For example, the Milky Way itself is a strong 
quadrupole source. 

Another measure of large scale structure in the Universe which depends explicitly on the power spectrum 
P(k) is the large scale bulk motion in the Universe (Lynden-Bell et al. 1988, Gorski 1991, Dekel et al. 1990). 
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A particularly useful quantity describing various models is the mean squared peculiar velocity (vp{R)), 
averaged over some scale R. 

In this paper, we derive explicit expressions for the peculiar velocity (wp(i?)) as well as the temperature 
correlation function C{0), both with and without its first few moments removed, as integrals over the mass 
autocorrelation function, ^(r) or its own integral J3(r). We also derive similar expressions for each term in 
the multipole expansion of C{6). In Section 2 we derive the basic relations, while in Section 3, we discuss 
how the various correlation functions interrelate. 



2 Fundamental Relations 
2.1 Definitions 

In any N-dimensional space, the correlation function is a useful statistic of quantifying some of the 
relationships between values of a function at different spatial positions. The function can be continuous, as 
in the ease of fluctuations in temperature or density, or it can even be discrete. A conceptually easy example 
of the latter is that of the autocovariance function of galaxies in three dimensions. Here it is a measure of 
the probability of finding a neighbor in the data set at a given distance from any galaxy (see Peebles, 1981, 
for example). 

The correlation function of any continuous distribution is defined mathematically in terms of an average 
over all of N-dimensional space. Hence, the angular correlation of temperature fluctuations in the CMBR 
on the scale 9 is defined as 

{ST\ST2) 

and the mass autocorrelation function for matter a distance r apart is given by 



C{6) = (1) 



= (2) 

The quantities 5T and 5p are respectively the temperature and density excess over the mean, with subscripts 
1 and 2 denoting position separated by an angle 9 or distance r. The angular brackets signify averaging in 
two dimensions over the sky and over all three dimensional space. It is assumed that the correlations are 
independent of direction. 

The power spectrum P{k) is a mcasin-e of the amplitude of density perturbations in Fourier space for 
a given wave- number. It is defined simply as \Sk\ , where 5k is the Fourier transform of the dimensionless 
density excess Sp/ {p). As a result of these definitions, it is obvious that P{k) is the three dimensional Fourier 
transform of $(r). Since the autocorrelation function is independent of direction, one can write 

P{k) = \6k\' =A^ drr^^rf-^. (3) 

The power spectrum (3) often arises in expressions for observable quantities which arise from density or 
potential fluctuations. Any such observable quantity, Y, usually involves integrals in Fourier space of the 
form 

/•oo 

(4) 



/•OO 

Y{x) = / dkk^ P{k)FY{k,x). 
Jo 



The quantity Fy (fc, x) explicitly relates the observed quantity Y (which may depend on some variable x) to 
the amplitude of waves having wavenumber k. 

We can relate the observed quantity Y to the autocorrelation function by substituting for P{k) from 

equation (3) and performing the fc-integral to yield an expression involving an integral in real space. In this 
way we derive the kernel, Gy(r, x), for mass fluctuations, peculiar velocities and temperature variations in 
the CMBR, such that 

Y{x)= I drr'^^{r)GY{r,x). (5) 
Jo 

This kernel is related directly to Fy by the expression 

Gy(r,a;) = — / dkk smikr) FyVk^x), (6) 
r Jo 
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or by the inverse expression 



1 

PY{k,x) = —^ drr sm{kr)GY{r,x). (7) 
Integration by parts of equation (5) leads to 

Y{x) = - dvMr)-^, (8) 

where ^_ 

Mr) = I drr^^ir). (9) 
Jo 

Since by definition, ^(r) must satisfy the integral eonstraint J:j(cc) = 0, while this is simply a boundary 
condition on J3, one can gain some insight into the underlying physics by considering (as in equation (8)) 
the derivative of the kernel dGy / dr rather than the kernel Gy itself. 

2.2 Mass Fluctuations and Peculiar Velocities. 

Consider a sphere having radius R. Then the mean squared mass excess in a randomly placed sphere 
of this radius is given by 

/fSM 



2 

V M 



(i?)M = ^^ dkk''P{k)W\kR), (10) 



where W{kR) is the window function, which is given by a Fourier transform of a three dimensional top hat 
sphere of radius R: 

^^'""^ = ^/pwe'^'"^'" = ^[sinfcE- .i^cosfcE] = 3^, (11) 
where ji is a spherical Bessel function. In analogy with eq. (5), we can rewrite equation (10) as 

{{w^^^y) = l^r'drar)GsM{r,R), (12) 
where equation (6) allows us to explicitly evaluate Gsm as 

GMr,R)-^,[^^ otherwise, ^^'^ 

where x = r/2R. Clearly, Gsm goes to zero for r > 2R since correlations larger than the diameter of the 
sphere are irrelevant. 

As an alternative to the top hat expression ((11)), one can smooth the power spectrum by weighting 
the points with a Gaussian window function of the form 

W{kRg) = e-'''</^, (14) 

where Rg is the Gaussian smoothing length. For this Gaussian window function, equation (6) allows us to 
obtain the kernel 

G.M(.i?)^^e-V..?. (15) 

In a sense, (15) is obvious since it is again the window function in real space. Some papers relate the Gaussian 
smoothing length to the top hat sphere radius. For example, Blumenthal et al. (1992) require the Gaussian 
window to have the same spatial volume as the top hat sphere: 

2 ^ 
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while Lccar (1993) points out that for Rg « A7R, the Gaussian window function gives similar integrals to 
the top hat for a wide range of power spectra. 

The kernels Gsm for both the top hat window function and for the Gaussian window with Rg given by 
(16) arc plotted in Figure la. Note that the expression for SM/M follows directly from the definition of £_{r) 
rather than from an evolutionary calculation. Thus the kernel for the top-hat window is zero for distances 
larger than the distance across the top-hat (2i?) since by definition, mass distributions are not sampled on 
these scales. The kernel for the Gaussian extends to cover all space since the weighted average is taken on 
all scales. In Figure lb we have plotted the derivative of the kernel, which is used with the J3 integral (8). 

In a similar fashion, one can calculate, to first order in perturbation theory, the rms peculiar velocity 
averaged over randomly placed regions of size R: 



fc2 



(17) 



where is well approximated by Q^'^ for density parameter < < 1 (Peebles 1976) and is approx- 

imately O"'/^ for « 1 (Lightman and Schcchter, 1991), and Q is the dimensionless density parameter. 
Actually, !F represents the logarithmic rate of change of the amplitude of fluctuations with scale factor. The 
expression (17) is valid when the cosmological constant in zero, and must be modified in the case of a 
nonzero cosmological constant. Ho is the current value of Bubble's parameter. For a top-hat window, (6) 
yields the peculiar velocity kernel: 

The comparable expression using the Gaussian window-function (14) is given by 

Gv(r, R) = ^o^^^^ erf (r/2Eg), (19) 

where erf is the error function. 

The velocity kernels Gv (and their derivatives) for both the top-hat window and the Gaussian window 
with i?g given by (16) are also plotted in Figure 1. Here, the difference between top hat and Gaussian is 
more pronounced in the kernel derivative (Figure lb) than in the kernel itself. This is a reflection of the fact 
that the autocorrelation function is in some sense a measure of the density distribution while J3 measures 
the mass distribution and it is on this that the peculiar velocity is more heavily dependent. The kernel 
falls as on large scales (r R). This behavior is obvious since the large scale contribution goes as 
{vp} oc M{r)/r oc pJ3{r)/r oc p(^{r)r^dr/r, and large regions can have a significant effect on the velocity, 
as can be seen already from equation (17). It is also similar to the expression for mean squared velocity in 
Peebles 1981 (equation 14.10), 

r), (20) 



{vl)=HlT\n) drr^{ 
Jo 



which explains the limiting behavior of Gv HqJ^'^(P)/2Rx. 

Kashlinsky (1992) used a similar method to relate both the peculiar velocity correlation function {y{r) = 
(vpj^.Vp2)) to ^(r) and the rms peculiar velocity to the projected angular correlation function of matter, w{6), 
found in the APM survey (Maddox et al. ,1990). 

2.3 Angular Correlation Function. 

The Sachs- Wolfe result for the CMBR angular correlation function is given by 







/ dhk"^ P{k)Fc{k,0) 
Jo 



(21) 
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where y = 2_R^sin(0/2) and Rh = 2c/ (Sachs and Wolfe, 1967). Strictly speaking, this is the correct 
expression only for the Einstein-deSitter f2 = 1, A = universe, where A is the cosmological constant. 
When A = and O 7^ 1, it is easy to generalize this for small enough angles corresponding to scales which 
crossed the horizon at early times when spatial curvature had little effect on the evolution of the Universe. 
Hence, for f2 < 1, expression (21) is a good approximation for small angles 9 < ^ radians. In most cases, 
the Sachs- Wolfe effect is the dominant contributor to the background fluctuations for 9 » 30' f2^/^; at these 
large angles, scattering effects are negligible. 

A slight variant of equation (21) also remains a good approximation in the case of a zero curvature 
model with A = 3(1 — Q.)H'^/c? (Kofman and Starobinsky 1985, Peebles 1984). Again, for small angles, 
Peebles (1984) showed that Rh should be replaced by the comoving distance 

'^'^w/"**"' + '''' 

and that the expression J-"(f2) should be replaced by an integral, which can be approximated as fJ^-^'' 
(Efstathiou et al. 1992). So that our results can be valid for all these cases, we will investigate the kernel 
Fc{k,9), defined in equation (21). 

Expanding Fc{k, 9) in terms of Legendre polynomials shows the magnitude of the multipole contribu- 
tions. 

00 00 
Fc{k,9) = A;-4^(2Z + 1) [ji{kRh)]^ Pi {cos 6) = J2 M^) Pi {cos 9), (23) 
1=0 1=0 
where ji are spherical bessel functions and Pi are Legendre polynomials. Figure 2 shows the amplitude of 
each of these contributions to the kernel (for fixed 6) as a function of wavenumber. The curves are each 
renormalized to a maximum of unity, and the monopole and dipole terms, which are singular at fc = 0, are 
not shown. Clearly the dominant contribution to the integral in (21) for small k (corresponding to large 
spatial scales) is from the monopole and dipole moments as all others are finite at the origin. As Figure 
2 shows, each subsequent, low order moment has a specific wavenumber associated with it, and picks out 
a small range of wavenumbers in the power spectrum P{k). In particular, the quadrupole contribution is 
dominated be the very shortest wavenumbers. The position of the maximum is a monotonically increasing 
function of while the amplitude decreases with /. The distance in fc-space between the maxima of adjacent 
/ tends to a constant value while the width of the peak increases. Hence the strong association of a given 
wavenumber with a given term in the expansion weakens significantly for higher order terms. 

In a similar fashion to equation (23), one can expand the spatial kernel in terms of Legendre polynomials: 

00 

Gc{e,T)^Y.9iir)Pi{e). (24) 

One normally ignores the monopole {I = 0) term since it represents only a renormalization of the mean 
temperature. Each of the gi{r) is related to the corresponding fi{k) from equation (23) using the integral 
(6). Defining x = r/2iih, the dipole term is given by 

^^(^^=^^{4, ^ ^ otherwise. ^''^ 

This term is normally ignored as well in cosmological comparisons since the major contributor to the observed 
dipole is thought to be the Doppler shift coming from the peculiar motion of our own galaxy. Doppler terms 
were explicitly excluded in Sachs and Wolfe's original calculation. The fact that contributions to the dipole 
can come from scales larger than the horizon lends credence to the argument (Paczynski and Piran, 1990) 
that it is possible for a significant part of the dipole moment not to be Doppler in origin. 

The higher order moments are all zero for a; > 1. For a; < 1 and Z > 1 we find the general result 



gi{r) = Rh 



1 1 2 

— ^a; 



2(2/-|-3)(2Z- 1) 6 



(2^ + 1) ' (-l)^(; + fc)! / 2 3_ 



12 ^ {l-k)\{k\Y \2k + l l + fc 2fc-|-3 k + 2 

fe=0 



2,2/s+3 



(26) 
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Figure 3a shows plots of these spatial kernels. The curves are again renormalized to a maximum of unity, 
and the monopole term, which is singular, is not plotted. All moments of higher order than the dipole are 
zero beyond the horizon and have {I — 1) zero points within it. As in the equivalent fc-space plot (figure 2) 
the maximum value decreases with increasing I. Here there is no specific distance associated with any given 
/ as all the peaks are at the origin. However, there is some suggestion to a characteristic range in r for each 
I as the width of the peaks decreases monotonically. 

Figure 3b plots the derivatives of the functions shown in figure 3a, which are relevant for the J3 integra- 
tion. Higher moments have maxima at progressively smaller scales, just as implied by the case in k space. 
Essentially, higher moments arise from correlations (or mass concentrations) at smaller spatial scales. Note 
that these curves have once again been renormalized and the amplitude of the initial peak actually decreases 
with I. Because correlations at a given angle 9 involve sums over many values of I, these plots (figures 3a 
and b) lead us to expect that such correlations do not arise from fluctuations on a specific associated spatial 
scale, but rather from similar ranges in r. 

One can also transform the general kernel, Fc{0, k) defined in equation (23). Defining s = sm9/2, the 
kernel for the correlation function with the monopole term removed (to get rid of singularities) is 



Gc{0,r) = Rh{ 



(1 — 2s )24j, otherwise. 



This agrees with a transform of the function derived from first principles by Traschen and Eardley (1986). 
The derivative function is given by 



(s- 1^)^-16' ifO<a;<s 

-1 + 2- 2! ifs<a;<l (28) 
- ( 1 - 2s^ ) , otherwise. 



Contour plots of Gc{6, r) and its derivative are shown in Figures 4a and 4b respectively, with the dipole 

term also removed. These figures demonstrate that there is very little tendency for correlations at a given 
angle to correspond to a given spatial distance r for either the integral involving ^(r) or the one involving 
Ji{r). Such a tendency may exist for very small angles of a few degrees or less, but for the most part, it is 
absent. 

This is further demonstrated in Figures 5 and 6. In Figures 5a and 5b we show slices through the contour 
maps in Figures 4a and 4b. Essentially all of the curves (for all angles shown) give significant contributions 
to C{0) from similar regions of space. Figures 6a and 6b present similar results but with the quadrupole 
term removed as well. Now the curves have more structure, but there is again little or no tendency for 
correlations on larger angles to arise from greater distances. 



3 Conclusions. 

We have derived explicit relations between the mass autocorrelation function ^(r) and the angular 
correlation of background radiation on the sky C{6) due to the Sachs- Wolfe effect. We also derived expressions 
relating C(0) (as well as its multipole moments) to the second moment of ^(r), namely Jsir). We argued that 
the integral involving Jz{r) is in some sense the more fundamental, since unlike the autocorrelation function, 
J3(r) does not satisfy an integral constraint. Because of the integral constraint, the kernel Gcid, r) does not 
provide the unique contribution from every radius r since every positive correlation must be balanced by a 
negative correlation somewhere else. Jair) does not suffer this problem. 

The expressions derived here may be used either with theoretical or observed mass autocorrelation 
functions to derive the microwave background fiuctuations. Of course, it is true that large galaxy surveys 
can be used to derive either an autocorrelation function or a power spectrum (Davis and Peebles 1983, 
Maddox et al. , 1990). Hence either equation (21) or equation (27) can in principle be used to derive the 
anisotropy. However, the spatial integral provides a direct measure of the volume that must be surveyed to 
obtain an accurate result at a given angle. If the dipole term (equation (25)) is also removed from the kernel 
it does not require integration over scales larger than the horizon. 
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The most surprising result that wc find is that mass fluctuations on all spatial scales contribute to 
temperature fluctuations on any but the very smallest angular scales, as is demonstrated in Figures 4 and 
5. This arises in part because the correlation at any angle is the sum of contributions from many multipole 
moments. Although there is a clear tendency for the higher multipole moments to arise from smaller distances 
r, when one adds the many multipole moments to obtain Gc, this tendency disappears. Hence it is dangerous 
and incorrect to assume that each angle 6 is related to a unique distance scale r for any correlation. 

We would like to thank Luiz DaCosta and Sandra Faber for helpful discussions. This work was supported 
in part by grants from NSF, NASA, and a UCSC faculty research grant. 
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Figure Captions 

Figure 1 a. Derived kernels for mass fluctuations and bulk velocities. Gsm plotted for a top-hat 
window function of radius R (solid line) and a Gaussian window of the same volume (dotted line). 
RGv/HgJ^^{n) plotted for a top-hat window (short dashed Hue) and a Gaussian window (long 
dashed line), b. Kernels for integration over J3 found by differentiating the functions shown in 
Figure la. 

Figure 2 Amplitude of multipoles of the fc-space kernel ,fi{k) for I = 2,3, ..6. These curves have been 
renormalized to a maximum of unity. 

Figure 3 a. Amplitude of multipoles for the real-space kernel, gi(r) for I = 1, 2.. 6. The curves have been 
renormalized to give a maximum of unity, b. Kernels for integration over J3 found by differentiating 
the functions shown in Figure 3a. The curves have been renormalized to give a maximum of unity. 

Figure 4 a. Contours of Gc{r,9)/Rf^ with monopole and dipole only subtracted. The maximum occurs 
at ^ = r = and has value j^. Negative contours are indicated by dashed lines. Levels are shown 
at logarithmic intervals with spacing of Adex with the zero level contour in bold and the lowest 
magnitude contour at 10^*'. b. Derivatives of function shown in Figure 4a. 

Figure 5 a. The value of the kernel Gc(^^. r)/i?)^ (with the monopole and dipole removed) for s = 
sin(6'/2) — 0.1, 0.3, 0.5, 0.7, 0.9 as functions of r. These are just slices through Figure 4a. b. Kernels 
for integration over J3 found by differentiating the functions shown in Figure 5a. 

Figure 6 a. The value of the kernel Gc(P,i")/R\ (with the monopole, dipole, and quadrupolc removed) 
for s = sin(0/2) = 0.1,0.3,0.5,0.7,0.9 as functions of r. This is the same as Figure 5a but with 
the quadrupole removed as well. b. Kernels for integration over J3 found by differentiating the 
functions shown in Figure 6a. 
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